Transcriptome organization of white blood cells through gene co-expression network analysis in a large RNA-seq dataset.

Gene co-expression network analysis enables identification of biologically meaningful clusters of co-regulated genes (modules) in an unsupervised manner. We present here the largest study conducted thus far of co-expression networks in white blood cells (WBC) based on RNA-seq data from 624 individuals. We identify 41 modules, 13 of them related to specific immune-related functions and cell types (e.g. neutrophils, B and T cells, NK cells, and plasmacytoid dendritic cells); we highlight biologically relevant lncRNAs for each annotated module of co-expressed genes. We further characterize with unprecedented resolution the modules in T cell sub-types, through the availability of 95 immune phenotypes obtained by flow cytometry in the same individuals. This study provides novel insights into the transcriptional architecture of human leukocytes, showing how network analysis can advance our understanding of coding and non-coding gene interactions in immune system cells.

Region-specific changes in gene expression are associated with cognitive deficits in the alpha-synuclein-induced model of Parkinson's disease: A transcriptomic profiling study.

Mild cognitive impairment (MCI) is a common trait of Parkinson's disease (PD), often associated with early motor deficits, eventually evolving to PD with dementia in later disease stages. The neuropathological substrate of MCI is poorly understood, which weakens the development and administration of proper therapies. In an α-synuclein (αSyn)-based model of PD featuring early motor and cognitive impairments, we investigated the transcriptome profile of brain regions involved in PD with cognitive deficits, via a transcriptomic analysis based on RNA sequencing (RNA-seq) technology. Rats infused in the substantia nigra with human α-synuclein oligomers (H-SynOs) developed mild cognitive deficits after three months, as measured by the two-trial recognition test in a Y-maze and the novel object recognition test. RNA-seq analysis showed that 17,436 genes were expressed in the anterior cingulate cortex (ACC) and 17,216 genes in the hippocampus (HC). In the ACC, 51 genes were differentially expressed between vehicle and H-αSynOs treated samples, which showed N= 21 upregulated and N = 30 downregulated genes. In the HC, 104 genes were differentially expressed, the majority of them not overlapping with DEGs in the ACC, with N = 41 upregulated and N = 63 downregulated in H-αSynOs-treated samples. The Gene Ontology (GO) and the Kyoto Encyclopedia of Gene and Genomes (KEGG) analysis, followed by the protein-protein interaction (PPI) network inspection of DEGs, revealed that in the ACC most enriched terms were related with immune functions, specifically with antigen processing/presentation via the major histocompatibility complex (MHC) class II and phagocytosis via CD68, supporting a role for dysregulated immune responses in early PD cognitive dysfunction. Immunofluorescence analysis confirmed the decreased expression of CD68 within microglial cells. In contrast, the most significantly enriched terms in the HC were mainly involved in mitochondrial homeostasis, potassium voltage-gated channel, cytoskeleton and fiber organisation, suggesting that the gene expression in the neuronal population was mostly affected in this region in early disease stages. Altogether results show that H-αSynOs trigger a region-specific dysregulation of gene expression in ACC and HC, providing a pathological substrate for MCI associated with early PD.

miR-365-3p mediates BCL11A and SOX6 erythroid-specific coregulation: A new player in HbF activation.

Hemoglobin switching is a complex biological process not yet fully elucidated. The mechanism regulating the suppression of fetal hemoglobin (HbF) expression is of particular interest because of the positive impact of HbF on the course of diseases such as ß-thalassemia and sickle cell disease, hereditary hemoglobin disorders that affect the health of countless individuals worldwide. Several transcription factors have been implicated in the control of HbF, of which BCL11A has emerged as a major player in HbF silencing. SOX6 has also been implicated in silencing HbF and is critical to the silencing of the mouse embryonic hemoglobins. BCL11A and SOX6 are co-expressed and physically interact in the erythroid compartment during differentiation. In this study, we observe that BCL11A knockout leads to post-transcriptional downregulation of SOX6 through activation of microRNA (miR)-365-3p. Downregulating SOX6 by transient ectopic expression of miR-365-3p or gene editing activates embryonic and fetal ß-like globin gene expression in erythroid cells. The synchronized expression of BCL11A and SOX6 is crucial for hemoglobin switching. In this study, we identified a BCL11A/miR-365-3p/SOX6 evolutionarily conserved pathway, providing insights into the regulation of the embryonic and fetal globin genes suggesting new targets for treating ß-hemoglobinopathies.

The Intranigral Infusion of Human-Alpha Synuclein Oligomers Induces a Cognitive Impairment in Rats Associated with Changes in Neuronal Firing and Neuroinflammation in the Anterior Cingulate Cortex.

Parkinson's disease (PD) is a complex pathology causing a plethora of non-motor symptoms besides classical motor impairments, including cognitive disturbances. Recent studies in the PD human brain have reported microgliosis in limbic and neocortical structures, suggesting a role for neuroinflammation in the development of cognitive decline. Yet, the mechanism underlying the cognitive pathology is under investigated, mainly for the lack of a valid preclinical neuropathological model reproducing the disease's motor and non-motor aspects. Here, we show that the bilateral intracerebral infusion of pre-formed human alpha synuclein oligomers (H-αSynOs) within the substantia nigra pars compacta (SNpc) offers a valid model for studying the cognitive symptoms of PD, which adds to the classical motor aspects previously described in the same model. Indeed, H-αSynOs-infused rats displayed memory deficits in the two-trial recognition task in a Y maze and the novel object recognition (NOR) test performed three months after the oligomer infusion. In the anterior cingulate cortex (ACC) of H-αSynOs-infused rats the in vivo electrophysiological activity was altered and the expression of the neuron-specific immediate early gene (IEG) Npas4 (Neuronal PAS domain protein 4) and the AMPA receptor subunit GluR1 were decreased. The histological analysis of the brain of cognitively impaired rats showed a neuroinflammatory response in cognition-related regions such as the ACC and discrete subareas of the hippocampus, in the absence of any evident neuronal loss, supporting a role of neuroinflammation in cognitive decline. We found an increased GFAP reactivity and the acquisition of a proinflammatory phenotype by microglia, as indicated by the increased levels of microglial Tumor Necrosis Factor alpha (TNF-α) as compared to vehicle-infused rats. Moreover, diffused deposits of phospho-alpha synuclein (p-αSyn) and Lewy neurite-like aggregates were found in the SNpc and striatum, suggesting the spreading of toxic protein within anatomically interconnected areas. Altogether, we present a neuropathological rat model of PD that is relevant for the study of cognitive dysfunction featuring the disease. The intranigral infusion of toxic oligomeric species of alpha-synuclein (α-Syn) induced spreading and neuroinflammation in distant cognition-relevant regions, which may drive the altered neuronal activity underlying cognitive deficits.

PRF1 mutation alters immune system activation, inflammation, and risk of autoimmunity.

BACKGROUND: Defective alleles within the PRF1 gene, encoding the pore-forming protein perforin, in combination with environmental factors, cause familial type 2 hemophagocytic lymphohistiocytosis (FHL2), a rare, severe autosomal recessive childhood disorder characterized by massive release of cytokines-cytokine storm. OBJECTIVE: The aim of this study was to determine the function of hypomorph PRF1:p.A91V g.72360387 G > A on multiple sclerosis (MS) and type 1 diabetes (T1D). METHODS: We cross-compare the association data for PRF1:p.A91V mutation derived from GWAS on adult MS and pediatric T1D in Sardinians. The novel association with T1D was replicated in metanalysis in 12,584 cases and 17,692 controls from Sardinia, the United Kingdom, and Scotland. To dissect this mutation function, we searched through the coincident association immunophenotypes in additional set of general population Sardinians. RESULTS: We report that PRF1:p.A91V, is associated with increase of lymphocyte levels, especially within the cytotoxic memory T-cells, at general population level with reduced interleukin 7 receptor expression on these cells. The minor allele increased risk of MS, in 2903 cases and 2880 controls from Sardinia p = 2.06 × 10-4, odds ratio OR = 1.29, replicating a previous finding, whereas it protects from T1D p = 1.04 × 10-5, OR = 0.82. CONCLUSION: Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

Author Correction: Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

We report on the influence of ~22 million variants on 731 immune cell traits in a cohort of 3,757 Sardinians. We detected 122 significant (P < 1.28 × 10-11) independent association signals for 459 cell traits at 70 loci (53 of them novel) identifying several molecules and mechanisms involved in cell regulation. Furthermore, 53 signals at 36 loci overlapped with previously reported disease-associated signals, predominantly for autoimmune disorders, highlighting intermediate phenotypes in pathogenesis. Collectively, our findings illustrate complex genetic regulation of immune cells with highly selective effects on autoimmune disease risk at the cell-subtype level. These results identify drug-targetable pathways informing the design of more specific treatments for autoimmune diseases.

Genetic regulation of gene expression and splicing during a 10-year period of human aging.

BACKGROUND: Molecular and cellular changes are intrinsic to aging and age-related diseases. Prior cross-sectional studies have investigated the combined effects of age and genetics on gene expression and alternative splicing; however, there has been no long-term, longitudinal characterization of these molecular changes, especially in older age. RESULTS: We perform RNA sequencing in whole blood from the same individuals at ages 70 and 80 to quantify how gene expression, alternative splicing, and their genetic regulation are altered during this 10-year period of advanced aging at a population and individual level. We observe that individuals are more similar to their own expression profiles later in life than profiles of other individuals their own age. We identify 1291 and 294 genes differentially expressed and alternatively spliced with age, as well as 529 genes with outlying individual trajectories. Further, we observe a strong correlation of genetic effects on expression and splicing between the two ages, with a small subset of tested genes showing a reduction in genetic associations with expression and splicing in older age. CONCLUSIONS: These findings demonstrate that, although the transcriptome and its genetic regulation is mostly stable late in life, a small subset of genes is dynamic and is characterized by a reduction in genetic regulation, most likely due to increasing environmental variance with age.

The proneural gene ASCL1 governs the transcriptional subgroup affiliation in glioblastoma stem cells by directly repressing the mesenchymal gene NDRG1.

Achaete-scute homolog 1 gene (ASCL1) is a gene classifier for the proneural (PN) transcriptional subgroup of glioblastoma (GBM) that has a relevant role in the neuronal-like differentiation of GBM cancer stem cells (CSCs) through the activation of a PN gene signature. Besides prototypical ASCL1 PN target genes, the molecular effectors mediating ASCL1 function in regulating GBM differentiation and, most relevantly, subgroup specification are currently unknown. Here we report that ASCL1 not only promotes the acquisition of a PN phenotype in CSCs by inducing a glial-to-neuronal lineage switch but also concomitantly represses mesenchymal (MES) features by directly downregulating the expression of N-Myc downstream-regulated gene 1 (NDRG1), which we propose as a novel gene classifier of MES GBMs. Increasing the expression of ASCL1 in PN CSCs results in suppression of self-renewal, promotion of differentiation and, most significantly, decrease in tumorigenesis, which is also reproduced by NDRG1 silencing. Conversely, both abrogation of ASCL1 expression in PN CSCs and enforcement of NDRG1 expression in either PN or MES CSCs induce proneural-to-mesenchymal transition (PMT) and enhanced mesenchymal features. Surprisingly, ASCL1 overexpression in MES CSCs increases malignant features and gives rise to a neuroendocrine-like secretory phenotype. Altogether, our results propose that the fine interplay between ASCL1 and its target NDRG1 might serve as potential subgroup-specific targetable vulnerability in GBM; enhancing ASCL1 expression in PN GBMs might reduce tumorigenesis, whereas repressing NDRG1 expression might be actionable to hamper the malignancy of GBM belonging to the MES subgroup.

Protective major histocompatibility complex allele prevents type 1 diabetes by shaping the intestinal microbiota early in ontogeny.

Certain MHC-II or HLA-D alleles dominantly protect from particular autoimmune diseases. For example, expression of the MHC-II Eα:Eß complex potently protects nonobese diabetic (NOD) mice, which normally lack this isotype, from spontaneous development of type 1 diabetes. However, the underlying mechanisms remain debated. We investigated MHC-II-mediated protection from type 1 diabetes using a previously reported NOD mouse line expressing an Eα transgene and, thereby, the Eα:Eß complex. Eα16/NOD females vertically protected their NOD offspring from diabetes and insulitis, an effect that was dependent on the intestinal microbiota; moreover, they developed autoimmunity when treated with certain antibiotics or raised in a germ-free environment. Genomic and proteomic analyses revealed NOD and Eα16/NOD mice to host mild but significant differences in the intestinal microbiotas during a critical early window of ontogeny, and transfer of cecal contents from the latter to the former suppressed insulitis. Thus, protection from autoimmunity afforded by particular MHC/HLA alleles can operate via intestinal microbes, highlighting potentially important societal implications of treating infants, or even just their pregnant mothers, with antibiotics.

Overexpression of the Cytokine BAFF and Autoimmunity Risk.

BACKGROUND: Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug-targetable pathways. METHODS: Using case-control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus-specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence-based fine mapping, cross-population and cross-phenotype analyses, and gene-expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. RESULTS: A variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease-risk allele was also associated with up-regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion-deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up-regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. CONCLUSIONS: A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.).

Population- and individual-specific regulatory variation in Sardinia.

Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic risk.

An Efficient Multiple-Testing Adjustment for eQTL Studies that Accounts for Linkage Disequilibrium between Variants.

Methods for multiple-testing correction in local expression quantitative trait locus (cis-eQTL) studies are a trade-off between statistical power and computational efficiency. Bonferroni correction, though computationally trivial, is overly conservative and fails to account for linkage disequilibrium between variants. Permutation-based methods are more powerful, though computationally far more intensive. We present an alternative correction method called eigenMT, which runs over 500 times faster than permutations and has adjusted p values that closely approximate empirical ones. To achieve this speed while also maintaining the accuracy of permutation-based methods, we estimate the effective number of independent variants tested for association with a particular gene, termed Meff, by using the eigenvalue decomposition of the genotype correlation matrix. We employ a regularized estimator of the correlation matrix to ensure Meff is robust and yields adjusted p values that closely approximate p values from permutations. Finally, using a common genotype matrix, we show that eigenMT can be applied with even greater efficiency to studies across tissues or conditions. Our method provides a simpler, more efficient approach to multiple-testing correction than existing methods and fits within existing pipelines for eQTL discovery.

Genome-wide association analyses based on whole-genome sequencing in Sardinia provide insights into regulation of hemoglobin levels.

We report genome-wide association study results for the levels of A1, A2 and fetal hemoglobins, analyzed for the first time concurrently. Integrating high-density array genotyping and whole-genome sequencing in a large general population cohort from Sardinia, we detected 23 associations at 10 loci. Five signals are due to variants at previously undetected loci: MPHOSPH9, PLTP-PCIF1, ZFPM1 (FOG1), NFIX and CCND3. Among the signals at known loci, ten are new lead variants and four are new independent signals. Half of all variants also showed pleiotropic associations with different hemoglobins, which further corroborated some of the detected associations and identified features of coordinated hemoglobin species production.

Transcriptome analysis reveals differential splicing events in IPF lung tissue.

Idiopathic pulmonary fibrosis (IPF) is a complex disease in which a multitude of proteins and networks are disrupted. Interrogation of the transcriptome through RNA sequencing (RNA-Seq) enables the determination of genes whose differential expression is most significant in IPF, as well as the detection of alternative splicing events which are not easily observed with traditional microarray experiments. We sequenced messenger RNA from 8 IPF lung samples and 7 healthy controls on an Illumina HiSeq 2000, and found evidence for substantial differential gene expression and differential splicing. 873 genes were differentially expressed in IPF (FDR<5%), and 440 unique genes had significant differential splicing events in at least one exonic region (FDR<5%). We used qPCR to validate the differential exon usage in the second and third most significant exonic regions, in the genes COL6A3 (RNA-Seq adjusted pval = 7.18e-10) and POSTN (RNA-Seq adjusted pval = 2.06e-09), which encode the extracellular matrix proteins collagen alpha-3(VI) and periostin. The increased gene-level expression of periostin has been associated with IPF and its clinical progression, but its differential splicing has not been studied in the context of this disease. Our results suggest that alternative splicing of these and other genes may be involved in the pathogenesis of IPF. We have developed an interactive web application which allows users to explore the results of our RNA-Seq experiment, as well as those of two previously published microarray experiments, and we hope that this will serve as a resource for future investigations of gene regulation in IPF.

Transcriptome analysis reveals differential splicing events in IPF lung tissue.

Idiopathic pulmonary fibrosis (IPF) is a complex disease in which a multitude of proteins and networks are disrupted. Interrogation of the transcriptome through RNA sequencing (RNA-Seq) enables the determination of genes whose differential expression is most significant in IPF, as well as the detection of alternative splicing events which are not easily observed with traditional microarray experiments. We sequenced messenger RNA from 8 IPF lung samples and 7 healthy controls on an Illumina HiSeq 2000, and found evidence for substantial differential gene expression and differential splicing. 873 genes were differentially expressed in IPF (FDR<5%), and 440 unique genes had significant differential splicing events in at least one exonic region (FDR<5%). We used qPCR to validate the differential exon usage in the second and third most significant exonic regions, in the genes COL6A3 (RNA-Seq adjusted pval = 7.18e-10) and POSTN (RNA-Seq adjusted pval = 2.06e-09), which encode the extracellular matrix proteins collagen alpha-3(VI) and periostin. The increased gene-level expression of periostin has been associated with IPF and its clinical progression, but its differential splicing has not been studied in the context of this disease. Our results suggest that alternative splicing of these and other genes may be involved in the pathogenesis of IPF. We have developed an interactive web application which allows users to explore the results of our RNA-Seq experiment, as well as those of two previously published microarray experiments, and we hope that this will serve as a resource for future investigations of gene regulation in IPF.

Gene signatures distinguish stage-specific prostate cancer stem cells isolated from transgenic adenocarcinoma of the mouse prostate lesions and predict the malignancy of human tumors.

The relevant social and economic impact of prostate adenocarcinoma, one of the leading causes of death in men, urges critical improvements in knowledge of the pathogenesis and cure of this disease. These can also be achieved by implementing in vitro and in vivo preclinical models by taking advantage of prostate cancer stem cells (PCSCs). The best-characterized mouse model of prostate cancer is the transgenic adenocarcinoma of the mouse prostate (TRAMP) model. TRAMP mice develop a progressive lesion called prostatic intraepithelial neoplasia that evolves into adenocarcinoma (AD) between 24 and 30 weeks of age. ADs often metastasize to lymph nodes, lung, bones, and kidneys. Eventually, approximately 5% of the mice develop an androgen-independent neuroendocrine adenocarcinoma. Here we report the establishment of long-term self-renewing PCSC lines from the different stages of TRAMP progression by application of the neurosphere assay. Stage-specific prostate cell lines were endowed with the critical features expected from malignant bona fide cancer stem cells, namely, self-renewal, multipotency, and tumorigenicity. Notably, transcriptome analysis of stage-specific PCSCs resulted in the generation of well-defined, meaningful gene signatures, which identify distinct stages of human tumor progression. As such, TRAMP-derived PCSCs represent a novel and valuable preclinical model for elucidating the pathogenetic mechanisms leading to prostate adenocarcinoma and for the identification of molecular mediators to be pursued as therapeutic targets.

Coexpression of CD49b and LAG-3 identifies human and mouse T regulatory type 1 cells.

CD4(+) type 1 T regulatory (Tr1) cells are induced in the periphery and have a pivotal role in promoting and maintaining tolerance. The absence of surface markers that uniquely identify Tr1 cells has limited their study and clinical applications. By gene expression profiling of human Tr1 cell clones, we identified the surface markers CD49b and lymphocyte activation gene 3 (LAG-3) as being stably and selectively coexpressed on mouse and human Tr1 cells. We showed the specificity of these markers in mouse models of intestinal inflammation and helminth infection and in the peripheral blood of healthy volunteers. The coexpression of CD49b and LAG-3 enables the isolation of highly suppressive human Tr1 cells from in vitro anergized cultures and allows the tracking of Tr1 cells in the peripheral blood of subjects who developed tolerance after allogeneic hematopoietic stem cell transplantation. The use of these markers makes it feasible to track Tr1 cells in vivo and purify Tr1 cells for cell therapy to induce or restore tolerance in subjects with immune-mediated diseases.

Lentiviral vector-based insertional mutagenesis identifies genes associated with liver cancer.

Transposons and γ-retroviruses have been efficiently used as insertional mutagens in different tissues to identify molecular culprits of cancer. However, these systems are characterized by recurring integrations that accumulate in tumor cells and that hamper the identification of early cancer-driving events among bystander and progression-related events. We developed an insertional mutagenesis platform based on lentiviral vectors (LVVs) by which we could efficiently induce hepatocellular carcinoma (HCC) in three different mouse models. By virtue of the LVV's replication-deficient nature and broad genome-wide integration pattern, LVV-based insertional mutagenesis allowed identification of four previously unknown liver cancer-associated genes from a limited number of integrations. We validated the oncogenic potential of all the identified genes in vivo, with different levels of penetrance. The newly identified genes are likely to play a role in human cancer because they are upregulated, amplified and/or deleted in human HCCs and can predict clinical outcomes of patients.